Spin-Orbital Separation in a Solvable Kugel-Khomskii Problem in One-Dimension 
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A Kugel-Khomskii model of spin- 1/2 spins and pseudospin-1/2 orbitals with nearest-neighbor 
interaction is rigorously shown to exhibit spin-orbital separation by means of a non-local unitary 
transformation. This transformation completely decouples the spins from the orbitals in such a way 
that the spins become paramagnetic while the orbitals form the solvable XXZ Heisenberg chain. 
The nature of various correlations is discussed. The significance of this model problem, in relation 
to the recently observed spin-orbital separation in Sr2Cu03, is noted. The spin-orbital separation 
for the orbital pseudospins greater than 1/2 is also discussed briefly. 
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Introduction. The multi-orbital Mott-Hubbard insula- 
tors are known to exhibit a rich variety of orbito- magnetic 
order due to the coupled nature of spin and orbital de- 
grees of freedom. The theoretical framework that one 
uses to discuss such problems goes by the name of Kugel- 
Khomskii models. These are a natural extension of the 
Heisenberg spin-exchange to orbitally degenerate cases 
where the spins and the orbitals get coupled via superex- 
change to form the spin-orbital models [1]. The phonons, 
through Jahn- Teller coupling, also participate in the sys- 
tems with orbital degeneracy. However, at the very least, 
one simply focuses on the electronic part, the leading 
term in energy, that gives rise to the Kugel-Khomskii 
models within second order perturbation theory in the 
limit of strong local repulsion (for effectively one parti- 
cle, electron or hole, per unit cell). 

There are a lot of materials, such as KCuFa, LaMnOa, 
MnFa, Sr2Cu03 (to name a few), where the Kugel- 
Khomskii models are directly applicable [1-3]. For ex- 
ample, in KCuFa, the Cu^"*" ions, having [Ar]3d^ config- 
uration in the octahedral crystal field of F~ ions, can be 
treated as having one hole in the twofold degenerate eg 
orbitals {dj.2_y2 and dj,2). It is a Mott insulator, and has 
been studied as a Kugel-Khomskii problem [4, 5]. There 
are also cases of a different type in which the effective 
problem looks like a spin-orbital model, although in ac- 
tual one may not be dealing with degenerate orbitals [6] . 
For instance, an odd-legged spin-1/2 tube can be de- 
scribed in terms of the effective spin-1/2 and chirality 
variables. Since chirality can be treated as a pseudospin- 
1/2 object, an effective model for such a spin-tube is a 
Kugel-Khmoskii model in one dimension [7-9]. Clearly, 
there is much interest in studying such spin-orbital mod- 
els with different motivations. 

In this Letter, we present a solvable Kugel-Khomskii 
model in one dimension (Id), whose most significant fea- 
ture is the 'spin-orbital separation', an effect similar to 
the spin-charge separation in the interacting Id electrons 
(Luttinger liquids). Interestingly, a recent report of spin- 
orbital separation in Sr2Cu03 [3] makes our (indepen- 
dent) theoretical findings experimentally relevant. 



In a Kugel-Khomskii model, the electron spin is de- 
scribed by the Pauli operators <?; , where I is the site index 
of a lattice. The operators for the orbitals depend upon 
the details thereof. Here, we consider twofold degenerate 
orbitals, in which case the operators in the orbital space 
can be treated as pseudospin-1/2 objects described by 
another set of Pauli operators, f^. For a system of spin- 
1/2 spins and pseudospin-1/2 orbitals, a typical Kugel- 
Khomskii problem has the following generic form. 

I 

+.13(^1 ■^i+i){ri,Ti+i}^,] (1) 

Here, the model Hamiltonian, Hj^j^, is written on a Id 
lattice, and it only has the nearest-neighbor exchange in- 
teractions. Of course, in general, the 11^^^^ could have 
further neighbor exchanges, and be living on any lat- 
tice. The exchange interactions, ^{1,2.3} i ^tre assumed 
to be antiferromagnetic, unless stated to be otherwise. 
The symbol {f/, t/^iIa denotes an anisotropic exchange, 
Arj^rj^j^-l-2 {t^Tj^-^ + Tj^^T;^) , between the orbitals. Ev- 
idently, the Hj^j^ is SU(2) symmetric for spins, and has 
U(l) symmetry for orbitals. Therefore, the total spin 
and the z-componcnt of total pseudospin are conserved. 
When A A' = 1, it becomes SU(2)xSU(2) symmetric. 
These are complex models to investigate theoretically, 
and have attracted a lot of attention. 

Of the many cases of Hj^j^, quite a few are special 
as they allow for exact analytical solution of some kind. 
Most notable is the SU(4) symmetric model correspond- 
ing to Ji = J2 = J3 and A = A' = 1 whose ground state 
energy and elementary excitations are known using Bethe 
ansatz [10, 11]. There are other interesting cases as well. 
For example, when Ji = J2 = SJa for A = A' = 1, the 
exact ground state is doubly degenerate with the spins 
and the orbitals separately forming pairwise singlets on 
alternate bonds. The same ground state also holds true 
for A = A' = and Ji = 2J3 and J2 = 3J3 [12-14]. The 
quantum phase diagram for general Ji / J3 and J2 / J3 has 
also been investigated using DMRG [15, 16]. Different 
variations of Hkk are also known to give some inter- 



2 



esting exactly solvable models. Notable ones are a sym- 
metric XY spin-orbital model which can be completely 
solved using Jordan- Wigner transformation [17], and an 
integrable dimerized chain with J3 = — 2 and {Ji, J2) al- 
ternately equal to (|, ^) or (5, |) on successive nearest- 
neighbor bonds [18]. We too have found a special case 
which is the subject of our discussion in this Letter. 

Model. We present the following spin-orbital Hamilto- 
nian in one dimension. We were led to this special model 
through an unexpected route starting with the infinite U 
Hubbard model [19, 20]. 
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Here, Jz and are the nearest-neighbor interaction pa- 
rameters. While Jz measures the strength of Ising inter- 
action between orbitals, J± couples the spin-exchange, 
Xi ;_|_-^ = {1 + ■ (Tj^j^)/2, with the orbital quantum fiuc- 
tuations (the XY part of the orbital-exchange) . It corre- 
sponds to the Hkk for Ji = A' = 0, J2 = J3 = J_l/4 
and A = 4J^/J_l. Evidently, the H is SU(2)xU(l) sym- 
metric, and conserves as well as the total spin, 
i (T;. Below we will show that the H of Eq. (2) is ex- 
actly solvable through a transformation which completely 
decouples spins from orbitals, such that the decoupled 
problems (of spins and orbitals) are separately solvable. 
Recently, there has been some experimental interest in 
the spin-orbital separation in quasi Id systems [3]. The 
model H presents a realistic and rigorous theoretical ex- 
ample of the same. 

An important property of the spin-exchange operators, 
Xi is that Xfi^^ = I, where I is the Identity opera- 
tor. It prompts us to view as a unimodular 'oper- 
ator phase' factor. Since it appears in conjunction with 
Tj^Tj^-^^ in _ff , we are led to think that one may suitably 
gauge-transform rj^ to rid the H of X^ operators. It 
turns out that we can indeed achieve this by the proce- 
dure invented in Ref. [19] by the present author, for the 
infinite-[/ Hubbard model in Id. Its applicability to the 
spin-orbital problem was also briefly noted there. The 
method exploits the identity Xf^^-^ = I, and assumes 
open boundary condition. It enables us to completely 
decouple the spins from orbitals by applying a suitably 
defined non-local unitary transformation on H. The de- 
tails are as follows. 

Consider the J± term in H (denoted below as H±), 



(3) 



1=1 



where Xi+i^i = (= X^ as it is Hermitian), and 

L is the total number of lattice sites. Here, we have ex- 
plicitly put in the open boundary condition in the lattice 
summation. Now, define on bond (1,2) a unitary opera- 
tor, 2 = P2 +P2X2.1, where = (1±t|)/2 are the 



projections operators for the orbitals states on site 2. By 
applying it on H±, we get 



= J±^{rir2 + T+T^) + (t+X^^^t^ + 

L-l 

1=3 
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(4) 



As a result of the transformation under lA-^ 2; three things 
have happened to H^. First, the spin-exchange oper- 
ator X^ 2 has vanished from the bond (1,2). Now we 

only have {JxT2 '^'''t'^'x^- Second, the operator, X23, 
on bond (2,3) has got transformed into a new operator 
"^2 3 " ^12^2 3- Thus, the Xy2 hasn't quite disap- 
peared from ifj^. Instead, it has been transferred to the 
bond (2,3). Unlike ^^23, the new operator ^"2 3 is not 
Hermitian. But it is unitary, which is important. Fur- 
thermore, we adopt the notation, X^2 — '^23- Therefore, 
X2 3-^3 2 ~ ^^3 2'^2 3 = And third, the interactions 
on the bonds beyond the bond (2,3) remain unaffected. 
These observations suggest that we may similarly trans- 
fer ^"2 3 to bond (3,4) and so on, and eventually get rid of 
all the spin-exchange operators in H± on an open chain. 

The strategy is to eat-up one-by-onc the spin-exchange 
factors on successive bonds. The unitary operator which 



does this can be defined as: U = Y[f=i^i where 



P, 



l+l 



(5) 



Here, X, 



Ll+l 



Xi 2X2 3 ' ' ' Xi is the string of ex- 
change operators, and Pj^-^ = (1 ±T;^;^)/2 are the orbital 
projectors. Clearly, the U is a highly non-local unitary 
operator. By transforming under U, we find that 



L-l 

1=1 



(6) 



This is a remarkable result with two notable features. 
One, the orbitals have been completely decoupled from 
the spins, and two, the spin-spin interaction has com- 
pletely vanished. This is a rigorous demonstration of the 
spin-orbital separation in a Kugel-Khomskii problem. 

The resultant orbital problem in Eq. (6) is an XY chain 
which is well-known to be exactly solvable, as it can be 
converted to a free fermion Hamiltonian using Jordan- 
Wigner transformation [21]. Moreover, the spins form 
an ideal paramagnet due to the lack of any spin-spin 
interaction in the transformed problem. This obviously 
makes H± a solved Kugel-Khomskii problem, where one 
exactly knows all the eigenvalues and eigenstates. Every 
eigenvalue of H± is exponentially degenerate (~ 2^) due 
to the free spins. Even the ground state of H± has an 
extensive entropy of L log 2. 
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We further note that W J^i ^ = Ez '^i and W+rf U 
= t[ . Since J^i '^i and operators are invariant under 
U, it enables us to write a more general model, solvable 
through the same spin-orbital decoupling procedure used 
for H± in Eq. (6). A simple and realistic modification 
that we do to H± is to add the nearest-neighbor orbital 
Ising term, X); '^f '^i+i- This gives us the model H of 
Eq. (2). Clearly, due to the invariance of t[ under U, the 
H also exhibits exact spin-orbital separation. That is, 

L-l 

1=1 

which is the XXZ Heisenberg problem, solvable by Bcthc 
ansatz [22-24]. Thus, the H is an exactly solvable Kugel- 
Khomskii problem. We can add more arbitrary terms 
of the type, F„({rf}) + K(Ei '^f , ^^f, E/ ^z"), to 
and still rigorously achieve spin-orbital separation us- 
ing lA. For example, we can certainly add a term like, 
— J2ii^'^i + V^Hj where h and rj are the fields for spins 
and orbitals, respectively. However, in general, the re- 
sulting problem may not always be analytically solvable. 

As a function of Jz (setting J± = 1), the exact ground 
state of H behaves as follows. For \ Jz \ < 1/2, the orbital 
part is critical with power law correlations and gapless ex- 
citations, akin to the XY case {Jz = 0). For Jz < —1/2, 
the ground state is ferro-orbital and the elementary or- 
biton excitations are gapped. Similarly, for Jz > 1/2, the 
orbital ground state is Neel ordered with gapped orbital 
excitations. In all these qualitatively different regions on 
Jz-line, the spins remain ideally paramagnetic. 

At this point, we also like to put our understanding 
of H± in perspective with the results in Ref. [15], where 
the Hkk for a = A' = has been investigated using 
DMRG. Particularly for Ji = 0, the point J2/J3 = 1 
(which is J2 — 1/4 in their notation) was identified as 
a transition point (see Fig. 1 in Ref. [15]), above which 
the ground state is a direct product of the ferromagnetic 
spins and the orbital fermi-sea (through Jordan- Wigner 
mapping), \F)s ® \JWfs)o- Here, F stands for ferro- 
magnetic, JWfs for the Jordan- Wigner fermi-sea, and 
the subscripts s and o indicate the spin and the orbital 
subsystems respectively. For J2/J3 < 1, they proposed 
a spin-dimerized antiferromagnetic (DAF) ground state, 
\DAF) s^\JWfs)o- Here, we make an important observa- 
tion that their transition point J2/ J3 = 1 is same as our 
H±. Therefore, J2/J3 = 1 is in fact a special point in the 
quantum phase diagram, hitherto unrealized, with exact 
solvability for the complete eigen-spectrum. The correct 
ground state at this special point is a highly entropic 
manifold of 2^ eigenstates, U{\si,S2, ■ ■ ■ , s l) ®\JW f s) 0} , 
where |si, S2, . . . , sl) denotes an Ising configuration of L 
spin-l/2s, with s; =t or | (for Z = 1 to L). 

Correlations. For the model H quite a few different 
correlation functions can be computed, or understood us- 
ing the results known for the decoupled problem. 



The simplest thing one can compute is the static spin 
susceptibility, Xs - iUiEi where /5 = 

l/k^T is the inverse temperature. Consider (J^i'^i) = 
Tr {p J2i ""f }i where p is the equilibrium thermal density 
matrix, e~^^ /Z, for H. For the decoupled Hamiltonian, 
WHU, the thermal density matrix is W pU = S?) p^, 
where I is the 2-^-dimensional identity matrix for the spin 
subsystem and p^ is thermal density matrix for the XXZ 
orbitals. Since erf is invariant under lA, the expec- 
tation (E/'''^) ~ ^^''s {E/ "'^1 = 0' as it ought to be 
for a paramagnet. Moreover, ((Ei*^^)^) = ^- Therefore, 
Xs ^ i/ksT. Here, tr^ denotes the trace over spins only. 
Similarly, tro for the orbitals, and Tr — tr^ trg- 

For the orbital pscudospins, the longitudinal suscepti- 
bility, xo - |{((E/ ^if) - (Li '^f )^}' can also be calcu- 
lated from the decoupled problem because the t[ opera- 
tors are invariant under U. The XXZ Heisenberg prob- 
lem has a very rich mathematical literature from which 
we can gladly quote the results for Xo- For the isotropic 
XXX orbital case, Xo ^ 1 + 2in(Tb/r) ^'-'^ small temper- 
atures [25]. For the anisotropic case, with 2Jz/J± > 1, 
the orbital excitations are gapped. Therefore, Xo will 
exhibit exponential behavior at low temperatures. How- 
ever, for 2Jz/J± < 1 when the excitations are gapless, 
the Bethe ansatz and field-theory treatment has shown 
that Xo ^ 7r(7r-9) sine at low temperatures. Here, 

6 = cos^^ {2Jz/J±), and c is some constant. The expo- 
nent ( — for < §, and equal to 2 otherwise [25]. 

Next we discuss the two-point orbital correlation func- 
tion, CQ^(r) = (rfT;^^) = {/5rf rf^^}. The invariance 
of operators under U implies C'^^(r) = tro {Pot[t[;^^} . 
The zz orbital correlation is thus exactly same as in the 
corresponding XXZ problem. For example, in the ground 
state of H±, it is Cg^{r) = — [-^] for odd integer val- 
ues of r and for even r [21]. For the full problem with 
non-zero Jz , it has been a terribly hard job to find amica- 
ble analytic forms of the correlation functions, but using 
field theoretic techniques, the asymptotic behavior has 
been predicted to be Co^(r) (— l)l''l-\/ln |r-|/|r| for the 
isotropic XXX case. Unlike the zz correlation function, 
the computation of xx and yy orbital correlations does 
not simplify because rj^ and rf' operators are are not- 
invariant under U. 

The presence of an external magnetic field, ^hj^i'^^i 
however, simplifies the computation of some ground state 
correlations by selecting the fully polarized spin state. 
For h > 0, the ground state wavefunction is \ijjg) = U\ tt 
. . . t)s ® \XXZ)o = I n • • ■ t). ® \XXZ), where | tt 
. . . t)s and \XXZ)o denote the fully polarized spin state 
and the XXZ orbital ground state, respectively. Notably 
the lA acts like an identity operator on the fully polarized 
spin state. The spin-spin correlation in this case is trivial, 
as all the spins are pointing in the same direction. One 
can also compute the xx and yy orbital correlations in 
some cases. For example, in the ground state of H±^ (with 
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h ^ 0), {i^,\TfTr^M = {i>,\TfTlM,) ^ (-1)1^1/^, 

as known for the XY model in the hmit of large r [26] . 
Surely, one can quote a lot more results for various calcu- 
lable objects, as the literature for the XXZ model is vast. 
But we stop it here, and briefly discuss the applicability 
of H, in a suitably adopted form, to Sr2Cu03 [3]. 

Remarks in relation to Sr2CuO'i. It is a ferro-orbital 
material with a large energy gap to orbiton excitations 
(about 2.5 eV from dx2_y2 to dxz)- This can be accounted 
for by a negative Jz, and a large crystal field, rj, that se- 
lects the ferro-orbital state. For the antiferromagnetic 
spins, one may include the nearest-neighbor antiferro- 
magnetic exchange. Thus, a minimal model for Sr2Cu03 
would be, H + Si ■ (T/+i + rj^^ t[' , with J > and 
iJzl,?? > J. The J and J± can be estimated from the 
spinon continuum width 0.6eV) and the orbiton ener- 
gies (bandwidth ~ .2eV). Although the J term presents a 
'formal' difficulty to rigorous spin-orbital separation (be- 
cause it is not invariant under U), the above model has 
in it the essential aspects of Sr2Cu03, and is worth in- 
vestigating. [From a mathematical standpoint, one could 
replace J^i ■ by ^if ^ or by Po T,i ' ^i+i^cT 
(where P~ = Yiii^ ~ '''f )/2 is a ferro-orbital projector), 
to rigorously obtain the spin-orbital decoupling under U. 
The ground states for these three cases would be similar. 
But obviously the first one is a realistic choice.] 

Orbital pseudo spin > 1/2. Finally, we discuss the spin- 
orbital separation in systems with more than two orbitals 
per site. The idea is to demonstrate that, in principle, 
this phenomenon can also occur when the orbital pseu- 
dospin > 1/2. As an example, consider the following gen- 
eralization of the H for an arbitrary orbital pscudospin, 
T: Ht^ EPA'TUi + J±Xu+,{f+'^f^-Y + h.c.)], 
where Tf and Tj^ are the angular momentum operators 
representing the orbital pseudospin on the Z*'' site, and 
fj^^T ^ (T;=^)2r^ -We find that Ht also exhibits com- 
plete spin-orbital decoupling under the unitary transfor- 
mation, Ut = Ylfj^i ^^t(.^^^ + 1)' where Urp{l,l -|- 1) = 
P;^"^-|-(l— F;^^) J ; Is 3 gencralization of Eq. (5). Here, 
Pf^ = ] — T)i{—T\i is the projector for the lowest eigen- 
state of ff . We can show that U^HtUj. = Y.i[J zT^fi^^ 
+ J±{f+'^'^f{^Y + '^■c-)]: similar to the H. Likewise, 
it is also valid for the more general forms of Ht, as Tf 
operators and di are invariant under Ut- 

Summary. We have found an exactly solvable one- 
dimensional Kugel-Khomskii model [H of Eq. (2)]. It 
exhibits spin-orbital separation, achieved by a unitary 
transformation (U) which systematically rids the H of the 
spin-exchange operators. The transformed problem has 
free spins and the XXZ model for the orbitals. While the 
obvious symmetry, SU(2)xU(l), of H implies the conser- 
vation of Si and rf only, in actual, all the spins 
that vanish from H under U are conserved. Combin- 



ing this with the integrability of the XXZ chain, one, in 
principle, has a complete knowledge of all the conserved 
operators of the H . This problem is argued to be rele- 
vant to the case of Sr2Cu03. We have also discussed the 
spin-orbital separation for the orbital pseudospins > 1/2, 
through a simple generalization of the H. 
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